Fractional loop delays in adaptive optics modeling and control.

This paper revisits the problem of optimal (minimum variance) control for adaptive optics (AO) systems when measurement and command applications are asynchronous, resulting in a non-integer servo loop delay. When not properly accounted for, such fractional delays may severely degrade the AO performance, especially in the presence of high-frequency vibrations. We present evidence of this performance degradation thanks to in-lab experimental measurements on the Gran Telescopio Canarias Adaptive Optics (GTCAO) system controlled with standard suboptimal linear quadratic Gaussian (LQG) controllers. A constructive, easy to implement LQG control design is then proposed and validated in a simulation for vibrations affecting the tip-tilt modes. Our methodology is very interesting because it allows a performance assessment for any linear controller in terms of variance, rejection transfer functions, power spectral densities, and stability margins. We also show how the continuous-time disturbance model can be derived from standard discrete-time disturbance data-based modeling.

Zonal-based high-performance control in adaptive optics systems with application to astronomy and satellite tracking.

This paper presents a model-based approach to adaptive optics (AO) control based on a zonal (i.e., pixelized) representation of the incoming atmospheric turbulence. Describing the turbulence on a zonal basis enables the encapsulation of the standard frozen-flow assumption into a control-oriented model. A multilayer zonal model is proposed for single-conjugate AO (SCAO) systems. It includes an edge compensation mechanism involving limited support, which results in a sparser model structure. To further reduce the computational complexity, new resultant zonal models localized in the telescope pupil are proposed, with AR1 or AR2 structures, that match the spatial and temporal cross-correlations of the incoming turbulence. The global performance of the resulting linear quadratic Gaussian (LQG) regulator is evaluated using end-to-end simulations and compared to several existing controllers for two different configurations: a very large telescope SCAO and low earth orbit satellite tracking. The results show the high potential of the new approach and highlight possible trade-offs between the performance and complexity.

Linear controller error budget assessment for classical adaptive optics systems.

Understanding limitations of adaptive optics (AO) systems is crucial when designing new systems. In particular, analyzing the potential of different controllers is of great interest for the upcoming AO systems of the very large telescopes (VLTs) and extremely large telescopes (ELTs). This paper thus details a complete error budget assessment formalism, based on analytic formulas involving the disturbance temporal power spectral density (PSD) and the controller transfer function, and is applicable to any linear controller. This formalism is presented here for the special case of classical AO systems, but can be extended to any closed- or open-loop, single- or multi-conjugated AO configuration. Special attention is paid to the "control-dependent" errors, the importance of which is directly related to the type of control used in the AO system. The proposed method is applied to a NAOS/VLT-type single conjugated AO system, using disturbance PSD derived from a simulated turbulence trajectory or estimated from wavefront sensor measurements, enabling the construction of detailed error budgets for an integrator and different linear quadratic Gaussian controllers. Application to ELT-sized systems is discussed in the conclusion.

Eye-pupil displacement and prediction: effects on residual wavefront in adaptive optics retinal imaging.

This paper studies the effect of pupil displacements on the best achievable performance of retinal imaging adaptive optics (AO) systems, using 52 trajectories of horizontal and vertical displacements sampled at 80 Hz by a pupil tracker (PT) device on 13 different subjects. This effect is quantified in the form of minimal root mean square (rms) of the residual phase affecting image formation, as a function of the delay between PT measurement and wavefront correction. It is shown that simple dynamic models identified from data can be used to predict horizontal and vertical pupil displacements with greater accuracy (in terms of average rms) over short-term time horizons. The potential impact of these improvements on residual wavefront rms is investigated. These results allow to quantify the part of disturbances corrected by retinal imaging systems that are caused by relative displacements of an otherwise fixed or slowy-varying subject-dependent aberration. They also suggest that prediction has a limited impact on wavefront rms and that taking into account PT measurements in real time improves the performance of AO retinal imaging systems.

First on-sky SCAO validation of full LQG control with vibration mitigation on the CANARY pathfinder.

Adaptive optics provides real time correction of wavefront disturbances on ground based telescopes. Optimizing control and performance is a key issue for ever more demanding instruments on ever larger telescopes affected not only by atmospheric turbulence, but also by vibrations, windshake and tracking errors. Linear Quadratic Gaussian control achieves optimal correction when provided with a temporal model of the disturbance. We present in this paper the first on-sky results of a Kalman filter based LQG control with vibration mitigation on the CANARY instrument at the Nasmyth platform of the 4.2-m William Herschel Telescope. The results demonstrate a clear improvement of performance for full LQG compared with standard integrator control, and assess the additional improvement brought by vibration filtering with a tip-tilt model identified from on-sky data, thus validating the strategy retained on the instrument SPHERE at the VLT.

Distributed Kalman filtering compared to Fourier domain preconditioned conjugate gradient for laser guide star tomography on extremely large telescopes.

This paper discusses the performance and cost of two computationally efficient Fourier-based tomographic wavefront reconstruction algorithms for wide-field laser guide star (LGS) adaptive optics (AO). The first algorithm is the iterative Fourier domain preconditioned conjugate gradient (FDPCG) algorithm developed by Yang et al. [Appl. Opt.45, 5281 (2006)], combined with pseudo-open-loop control (POLC). FDPCG's computational cost is proportional to N log(N), where N denotes the dimensionality of the tomography problem. The second algorithm is the distributed Kalman filter (DKF) developed by Massioni et al. [J. Opt. Soc. Am. A28, 2298 (2011)], which is a noniterative spatially invariant controller. When implemented in the Fourier domain, DKF's cost is also proportional to N log(N). Both algorithms are capable of estimating spatial frequency components of the residual phase beyond the wavefront sensor (WFS) cutoff frequency thanks to regularization, thereby reducing WFS spatial aliasing at the expense of more computations. We present performance and cost analyses for the LGS multiconjugate AO system under design for the Thirty Meter Telescope, as well as DKF's sensitivity to uncertainties in wind profile prior information. We found that, provided the wind profile is known to better than 10% wind speed accuracy and 20 deg wind direction accuracy, DKF, despite its spatial invariance assumptions, delivers a significantly reduced wavefront error compared to the static FDPCG minimum variance estimator combined with POLC. Due to its nonsequential nature and high degree of parallelism, DKF is particularly well suited for real-time implementation on inexpensive off-the-shelf graphics processing units.

Fast computation of an optimal controller for large-scale adaptive optics.

The linear quadratic Gaussian regulator provides the minimum-variance control solution for a linear time-invariant system. For adaptive optics (AO) applications, under the hypothesis of a deformable mirror with instantaneous response, such a controller boils down to a minimum-variance phase estimator (a Kalman filter) and a projection onto the mirror space. The Kalman filter gain can be computed by solving an algebraic Riccati matrix equation, whose computational complexity grows very quickly with the size of the telescope aperture. This "curse of dimensionality" makes the standard solvers for Riccati equations very slow in the case of extremely large telescopes. In this article, we propose a way of computing the Kalman gain for AO systems by means of an approximation that considers the turbulence phase screen as the cropped version of an infinite-size screen. We demonstrate the advantages of the methods for both off- and on-line computational time, and we evaluate its performance for classical AO as well as for wide-field tomographic AO with multiple natural guide stars. Simulation results are reported.

Minimum-variance control for woofer-tweeter systems in adaptive optics.

The woofer-tweeter concept in adaptive optics consists in correcting for the turbulent wavefront disturbance with a combination of two deformable mirrors (DMs). The woofer corrects for temporally slow-evolving, spatially low-frequency, large-amplitude disturbances, whereas the tweeter is generally its complement, i.e., corrects for faster higher-order modes with lower amplitude. A special feature is that in general both are able to engender a common correction space. In this contribution a minimum-variance solution for the double stage woofer-tweeter concept in adaptive optics systems is addressed using a linear-quadratic-Gaussian approach. An analytical model is built upon previous developments on a single DM with temporal dynamics that accommodates a double-stage woofer-tweeter DM. Monte Carlo simulations are run for a system featuring an 8×8 actuator DM (considered infinitely fast), mounted on a steering tip/tilt platform (considered slow). Results show that it is essential to take into account temporal dynamics on the estimation step. Besides, unlike the other control strategies considered, the optimal solution is always stable.

Wide field adaptive optics laboratory demonstration with closed-loop tomographic control.

HOMER, the new bench developed at ONERA devoted to wide field adaptive optics (WFAO) laboratory research, has allowed the first experimental validations of multi-conjugate adaptive optics (MCAO) and laser tomography adaptive optics (LTAO) concepts with a linear quadratic Gaussian (LQG) control approach. Results obtained in LTAO in closed loop show the significant gain in performance brought by LQG control, which allows tomographic reconstruction. We present a calibration and model identification strategy. Experimental results are shown to be consistent with end-to-end simulations. These results are very encouraging and demonstrate robustness of performance with respect to inevitable experimental uncertainties. They represent a first step for the study of very large telescope (VLT) and extremely large telescopes (ELT) instruments.

On the optimal reconstruction and control of adaptive optical systems with mirror dynamics.

In adaptive optics (AO) the deformable mirror (DM) dynamics are usually neglected because, in general, the DM can be considered infinitely fast. Such assumption may no longer apply for the upcoming Extremely Large Telescopes (ELTs) with DM that are several meters in diameter with slow and/or resonant responses. For such systems an important challenge is to design an optimal regulator minimizing the variance of the residual phase. In this contribution, the general optimal minimum-variance (MV) solution to the full dynamical reconstruction and control problem of AO systems (AOSs) is established. It can be looked upon as the parent solution from which simpler (used hitherto) suboptimal solutions can be derived as special cases. These include either partial DM-dynamics-free solutions or solutions derived from the static minimum-variance reconstruction (where both atmospheric disturbance and DM dynamics are neglected altogether). Based on a continuous stochastic model of the disturbance, a state-space approach is developed that yields a fully optimal MV solution in the form of a discrete-time linear-quadratic-Gaussian (LQG) regulator design. From this LQG standpoint, the control-oriented state-space model allows one to (1) derive the optimal state-feedback linear regulator and (2) evaluate the performance of both the optimal and the sub-optimal solutions. Performance results are given for weakly damped second-order oscillatory DMs with large-amplitude resonant responses, in conditions representative of an ELT AO system. The highly energetic optical disturbance caused on the tip/tilt (TT) modes by the wind buffeting is considered. Results show that resonant responses are correctly handled with the MV regulator developed here. The use of sub-optimal regulators results in prohibitive performance losses in terms of residual variance; in addition, the closed-loop system may become unstable for resonant frequencies in the range of interest.

Linear quadratic Gaussian control for adaptive optics and multiconjugate adaptive optics: experimental and numerical analysis.

We present a comprehensive analysis of the linear quadratic Gaussian control approach applied to adaptive optics (AO) and multiconjugated AO (MCAO) based on numerical and experimental validations. The structure of the control law is presented and its main properties discussed. We then propose an extended experimental validation of this control law in AO and a simplified MCAO configuration. Performance is compared with end-to-end numerical simulations. Sensitivity of the performance regarding tuning parameters is tested. Finally, extension to full MCAO and laser tomographic AO (LTAO) through numerical simulation is presented and analyzed.

First laboratory validation of vibration filtering with LQG control law for adaptive optics.

We present a first experimental validation of vibration filtering with a Linear Quadratic Gaussian (LQG) control law in Adaptive Optics (AO). A quasi-pure mechanical vibration is generated on a classic AO bench and filtered by the control law, leading to an improvement of the Strehl Ratio and image stability. Vibration filtering may be applied to any AO system, but these results are of particular interest for eXtrem AO, and for instance for the SPHERE AO design, where high performance is required.

Optimal control, observers and integrators in adaptive optics.

The fundamental issue of residual phase variance minimization in adaptive optics (AO) loops is addressed here from a control engineering perspective. This problem, when suitably modeled using a state-space approach, can be broken down into an optimal deterministic control problem and an optimal estimation problem, the solution of which are a linear quadratic (LQ) control and a Kalman filter. This approach provides a convenient framework for analyzing existing AO controllers, which are shown to contain an implicit phase turbulent model. In particular, standard integrator-based AO controllers assume a constant turbulent phase, which renders them prone to the notorious wind-up effect.

Optimal control law for classical and multiconjugate adaptive optics.

Classical adaptive optics (AO) is now a widespread technique for high-resolution imaging with astronomical ground-based telescopes. It generally uses simple and efficient control algorithms. Multiconjugate adaptive optics (MCAO) is a more recent and very promising technique that should extend the corrected field of view. This technique has not yet been experimentally validated, but simulations already show its high potential. The importance for MCAO of an optimal reconstruction using turbulence spatial statistics has already been demonstrated through open-loop simulations. We propose an optimal closed-loop control law that accounts for both spatial and temporal statistics. The prior information on the turbulence, as well as on the wave-front sensing noise, is expressed in a state-space model. The optimal phase estimation is then given by a Kalman filter. The equations describing the system are given and the underlying assumptions explained. The control law is then derived. The gain brought by this approach is demonstrated through MCAO numerical simulations representative of astronomical observation on a 8-m-class telescope in the near infrared. We also discuss the application of this control approach to classical AO. Even in classical AO, the technique could be relevant especially for future extreme AO systems.